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SUMMARY 
STATEMENT OF PROBLEM 


Determine some of the theoretical properties of the multipole sound fields generated 
by both free turbulence and turbulence in the presence of solid boundaries at distances that 
are less than a wavelength from the turbulence. Compare the frequency spectra and spatial 
properties of the induction near-sound fields with the frequency spectra and spatial properties 
of the far-sound fields of these multipole sound sources. 


FINDINGS 


There are significant differences in the frequency spectra and the spatial properties 
of the various multipole components of the sound field at distances greater than and less than 
a wavelength from a turbulent region. Approximate corrections to datataken in the near-sound 
field of multipole sound sources can be made if the geometric extent, type, and orientation of 
the multipole sources are known or can be estimated. 


RECOMMENDATIONS 


If possible, measurements of the sound from multipole sound sources should be taken 
at a distance that is many wavelengths and many ‘‘diameters’’ from the sound sources. If 
this is not possible or desirable, information about the geometric extent, type, and orientation 
of the sound sources should be obtained. 
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NOTATION 


Velocity of sound in medium outside turbulent region 
Drag coefficient 


Instantaneous applied force or acoustic dipole strength per 
unit volume at point 7 in the turbulent region 

F.=F,G, t-r/c) =f 
f; iho (7, t-1r/c) 


Perpendicular distance from point to plane 


i 


Acoustic intensity in the 2, coordinate direction at point 
Acoustic intensity in the r, 9, and ¢ directions at point z 
Mass rate density of sources 

Typical linear dimension of the flow 


Mach number based on c and the convection velocity of the 
forces or the stresses 


Vector Mach number based on ¢ and the convection velocity of the 
forces or the stresses 


ith component of the vector Mach number 

M, = U/c, typical Mach number of the flow 

Sound power radiated by the turbulent region 
Instantaneous sound pressure at point & 

Mean-square sound pressure at point 7 

Reynolds number 

r = |z — y|, spherical polar coordinate centered at point 7 
7, = &; — y,, ith component of 7 

THe-y 
Correlation area 


Element of area at 7 


Instantaneous applied stress or acoustic quadrupole strength 
per unit volume at point 7 in the turbulent region 


T;, = T;; (¥, t-1/c) = tii 
hang BS OF (y, t—r/c) 
Time 


Typical velocity in the flow 


| 


2, j,k, 1,m 


u=u-7?/r, the component of @ in the 7-direction 

Vector velocity at 7 at retarded time ¢-r/c 

u, =u, (y, t-r/c), fluid velocity at 7 at retarded time ¢-r/c 
Correlation volume 

kth component of the instantaneous particle velocity at % 
Vector position of a point outside of the turbulent region 
ith coordinate of a point outside of the turbulent region 
Vector position of a point in the turbulent region 

ith coordinate of a point in the turbulent region 


Element of volume at 7 


Unit diagonal tensor 
Efficiency of conversion of mechanical to acoustic energy for dipole 


Efficiency of conversion of mechanical to acoustic energy for 
quadrupole 


na yt ie 
6 = |6|, spherical polar coordinate 

Vector in the direction of increasing 0 

ith component of 6 

Mean density of the medium 

Instantaneous velocity potential at z 

¢ = ||, spherical polar coordinate 

Vector in the direction of increasing ¢ 

ith component & 

Angular frequency or 27 times the frequency of the radiated sound 


Angular frequency of the velocity fluctuations 


Subscripts which take on values 1, 2, or 3 for the three coordinate 
directions 


Primes indicate condition at independent position 7’ in the 
turbulent region 


Dots indicate partial derivatives with respect to time 
Overbar on individual quantities indicates a vector 


Overbar on products of time-dependent quantities indicates a 
time average 


ABSTRACT 


General expressions for the velocity potential, particle velocity, instan- 
taneous sound pressure, mean-square sound pressure, and sound intensity for 
the sound generated by turbulence are derived, starting from Lighthill’s funda- 
mental relations for the density fluctuations in the medium outside a turbulent 
region. The dipole radiation from the turbulent boundary layer of a rigid bound- 
ary, the quadrupole radiation from isotropic turbulence, and the lateral quadrupole 
radiation from turbulence in the presence of a large mean shear are discussed. 
The frequency spectra and directional patterns of the sound observed at points 
in the sound field that are much less than a wavelength from multipole sound 
sources are shown to differ significantly from the frequency spectra and direc- 
tional patterns of the sound observed at points in the sound field that are many 


wavelengths from the multipole sound sources. 


INTRODUCTION 


Significant progress has been made toward an understanding of the subject of turbulence 
noise since the publication of Lighthill’s fundamental relations! for the density fluctuations 
in the medium outside a turbulent region. The expressions defining the radiated or far-sound 
field of a turbulent region in terms of multipole fields have been extensively discussed by 
Lighthill,? Proudman,® Phillips,* and others.°-® The purpose of this report is to present 
expressions which also explicitly define the near-sound field at distances either not large 
compared with a wavelength of the radiated sound or not large compared with the geometric 
extent of the turbulent region. It is hoped that these expressions will be useful in interpreting 
noise data taken close to turbulent jets and boundary layers. 

Acoustic multipoles are also associated with phenomena other than turbulence. The 
general expressions for the sound fields of multipoles, of course, are also applicable to these 
sound fields. For instance, the near-sound fields produced by rotating propellers® and by 
surface disturbances!° have been investigated and analyzed using acoustic dipole theory. 

Because the far-sound field is simpler and usually of more interest, measurements in 
the near-sound field of a turbulent region have usually been avoided. Some reports on the 
sound from turbulent jets, however, do contain data on the sound pressures in the near 
field.11—18 In general, the data indicate an enhancement of the low-frequency portion of the 
sound spectra and a different sound pressure-distance dependence than has been observed 
for the far-sound field. In the future, careful measurements of the noise close to turbulent 
boundary layers may alsoshow near-field effects.1? Since observations in the near-sound 
field are often necessary or convenient, it is advantageous to have means of interpreting the 
data. This report is intended to facilitate the making of corrections to existing and future 
data and the planning of new experiments. 


Dee ferences are listed on page 33. 


GENERAL SOUND FIELDS OF DIPOLES AND QUADRUPOLES 
DEFINITION OF NEAR FIELD 
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As used in this report, the unqualified expression ‘‘near-sound field’’ refers either to 
the region called the ‘‘induction near field’’ at distances less than a wavelength from the 
multipole source, or to the region called the ‘‘geometric near field’’ at distances from the 
sound source that are less than the order of the geometric extent of the source. In the induc- 
tion near field the particle velocity and pressure fields are not in phase in time. These fields 
contain a relatively large amount of unradiated energy that alternates between the kinetic 

and potential forms in the immediate vicinity of the multipole source. The geometric near 
field exists because an average distance from the measuring point to the source is usually 
used in presenting data for the region close to an extended source, even though the actual 
distances to the various parts of the source should be used. 

In many practical situations the regions defined by the two types of near field overlap 
and are comparable in extent, so that the effects of being in each region separately may be 
difficult to determine. The situation is further complicated by the fact that the effects of 
being in each region may tend to compensate each other. For example, if the distance from 
the nearest part of an extended multipole source is used as the characteristic distance from 
the source, 1! the rms sound pressure in the induction near field varies inversely as the square 
or some higher power of the distance from the source, and in the geometric near field the rms 


sound pressure may be nearly independent of the distance from the source. 


SOUND FIELD OUTSIDE A FIELD OF FLUCTUATING 
FORCES AND STRESSES 


The general expressions derived by Lighthill for density fluctuations in the medium 
outside a turbulent region have been expanded and applied, up to now, only for the region many 
wavelengths from the turbulent region. The purpose of this report is to expand these general 
expressions also for the region defined by the induction near field. As has been done for the 
radiation field, a number of special cases will also be discussed because under some practical 
conditions one type of multipole may predominate sufficiently to permit an approximate descrip- 
tion of the sound field by one set of simplified equations. 

By using the relation between density fluctuations and pressure fluctuations, the instan- 
taneous sound pressure p at time ¢ and point z outside aregion of fluctuating forces and a re- 
gion of fluctuating stresses, as sketched in Figure 1, is given, respectively, by the volume 


integrals 


1 a {F; a : 
p = - — |—|— ]ady (Equation 11, Reference 1) [1] 
4a Jou,\r 


Figure 1 — Location of Origin and Volume 
Element in a Turbulent Region 


and 


1 02 ee : 
== — } dz (Equation 15, Reference 1) [2] 
4a J Og, Oz. r 


where F’(¥, 2) is the zth component of the instantaneous applied force or acoustic dipole 
strength per unit volume at time ¢ and at point 7 in the turbulent region, 


T;; (y, ¢) is the 27 component of the instantaneous applied stress or acoustic quadrupole 
strength per unit volume at time ¢ and at point y in the turbulent region, 


F. = F.(y,t-1/c), 

T;; = T;; (y, t-7r/c), * 

ip = |z-y|, and 

Cc is the velocity of sound in the medium. 


If the applied forces are distributed over a surface, * Fis considered to be a surface density 
and the integral is a surface integral instead of a volume integral. 

If the differentiation in Equation [1] is performed, the general expression for the instan- 
taneous sound pressure due to the distribution of dipoles becomes 


ioe Be ay [3] 
4nJ 7 \ re 2 


where r, = a, -—y, and dots indicate partial differentiation with respect to time. The term 

in F; represents the radiated sound pressure, and the term in F, represents induction near-field 
sound pressure. The effect of the geometric near field is manifested in the factor 7,/r and in 
the dependence on r. In other words, if 7, can be approximated by z,, geometrical near-field 


effects are negligible. 


*This is a slight variation of the conventional notation used by Lighthill, necessitated by the need for a short 
symbol. 


Since the sound intensity, defined as the time average of the product of the particle 
velocity and the sound pressure at a point, cannot be approximated in the induction near field 
by D*/po where p, is the mean density and the overbar on the squared quantity signified a 
time average, it is necessary to obtain an expression for the particle velocity. This is most 
easily obtained by taking the gradient of the velocity potential ® which is first obtained by 


integrating the equation 


° 


P=- py? [4] 


After neglecting the constants of integration, the velocity potential for the distribution of 


dipoles is 


1 aff Pa Tg \ 
® =- |e —+—/]dy [5] 
AmpoJj 7 \re 2 


. 


where f, = f; (Y, t—r/c) and f, = F,.* The kth component of the particle velocity v or the 


particle velocity in the z, coordinate direction is then 


1 ope Wy OU | Bile Me i 
On = J Ss | Oe || Se OT [6] 
47 Po r2 \pre? re 73 ro 


where 65,, is the unit diagonal tensor and repeated subscripts are to be summed (z = 1, 2, 3). 


The &th component of the vector sound intensity is then given by 


Ppcrcgil OR NRO a Biaeaiy 4 

1 ij i i j j j 

ee eee et 
16 npg rr \TC 92 I\p%e2 p23 


He Pie te 
= — .)( + ) dy dy’ (7] 
r PO PY 08 Be 


where primes and different subscripts are used in the expression for the particle velocity to 


indicate functions of 7’. Similarly, the mean-square sound pressure is given by 


ser 1 re F. F. Fe Fe 
aoe | euy t say in [8] 
2 , TC 2 , 22 


rT r PRO -— 


1 

* 

-_-—_ ff fi ¥y, t)dy= B; (t), the dipole strength as usually defined. 
™Po 


4 


In like manner, if the double differentiation in Equation [2] is performed, the general 


expression for the instantaneous sound pressure due to the distribution of quadrupoles becomes 


Meese Mie 8 Tos 8 Te Fea 
p=— eat ee Leys = Oa: i pea dy [9] 
47 r? re? re re INr2g 73 


The term in ih represents the radiated sound pressure, and the terms in T represent the in- 
duction near-field sound pressure. The terms in T;; are Significant only in the transition 
region between the induction near field and the far field, and are conveniently lumped with 


the induction near-field terms. After again neglecting the constants of integration, the veloc- 


ity potential is 
i BVT Se WONT nee S) 6H os OND 
Ga [ep oe -~8,,(—-+—)la7 10 
CPD! || A 5 2?) 2 3 LT \\ 9) 3 


c = @ r TChin: 


where ti = ti (y, t-r/c) and t;, = T;,- This velocity potential, as well as the previous one, 
satisfies the wave equation 


; eryies 
vya=-— © {11] 
pe 


The &th component of the particle velocity is 


1 fs Sih, aD, a 
1) + + — 


kom ae 3 
47pg r3 rc? ro? rc 74 
r PANG T Tanne Bee 
k 
“ Ce 22, | ee | dy [19] 
r TEN 2) 2 noe 74 


By again using primes and different subscripts in the expression for the particle velocity and 


taking the time average of the product of p and v,, the kth component of the sound intensity 
becomes 


ee 


’ ee e eo 


re 2 Ea f i 
} 1 | iil ae ee ts) ( TAOS SET coca! 018 is 
k = a + — + + 


= ——]\ — + + 

16 7p 273 \pG2 r¢ 7 r’e3 2 Q2 r’73¢ 74 
rT r Ben AGE 37. 3 T. T, Bi) Bae 
tj m ij i i lm lm l 

- (0, = +2, 7\(4. —\( + + =| [13] 
r pe r° I \ re? 6 73 r°2g2 9736 mia 


Tatton T;; T;; ie 67," 15 ih 1507 has 
: + dy dy [13] 
roo p42 G2 r’3e pees 


Similarly, the mean-square sound pressure is given by 


= 1 fj rately Ui ( He) 8 ee 88 =a) i) BIE 8 = 
i = ———— | — + + + + 
16 7? 72’? 2 7? 3 ‘ce? 02 £3 


rC Cc rc r’“o r 
Gee Tj T;; T, ‘ 8 Th 3 Tim 
= 85, — + — — + 
2 ro or? ee p? 46 73 
Bin fTon  BTHe CORN Te aT 
ij ij ij ij m m 
2G 2 ee 
7? rez re r pIa@ Gee 
Te TTC 
l l 
+8558 1m i a ( — + =| dy dy” [14] 
me r’2¢ 3 


SIMPLIFYING ASSUMPTIONS 


These are indeed formidable expressions. A few reasonable assumptions may be made, 
however, to reduce themto more tractable forms. First, it may be assumed that the turbulence 
is sufficiently steady over intervals of time required for sound to travel over distances for 
which the velocity fluctuations are appreciably correlated so that simultaneous correlations 
may be taken at 7 and 7’. This means that the wavelength of the sound in question must be 
long compared to the ‘‘eddy size.’? Observed average eddy sizes are approximately a wave- 
length for frequencies of the order of ten times the frequency of the peak in observed noise 
spectra of turbulent jets. Since it is reasonable to assume that smaller than average eddies 
contribute most to the high-frequency components, the assumption of simultaneous correlations 
seems to be justified. If, in addition, it is assumed that ris much greater than the eddy diam- 
eter, then r= 7. 

Another reasonable assumption that permits a considerable reduction of the equations 
is that the fluctuating forces and fluctuating stresses are approximately stationary random 
functions of time. In other words, the turbulence is assumed to be essentially steady and not 
decaying appreciably.* This assumption still permits the existence of both gradual space 


so proudmanc has shown that, for the case of isotropic turbulence, the ‘‘decay’’ terms are indeed negligible in 


comparison with the ‘‘instantaneous’”’ terms. 


variations and gradual time variations in the mean values of the turbulent fluctuations. Since 


the operations of differentiating and averaging permute, this assumption is expressed by the 


equation??? 


me ene [g(@)]=0 [15] 
Bi te Ene i 


Therefore, reductions of the following type may be made: 


Ti, es iP, a= [16] 
and 
T;; qT), =0 [17] 


SIMPLIFIED GENERAL EXPRESSIONS FOR THE SOUND FIELD 


By making the above assumptions, Equations [7], [8], [13], and [14] become 


1 "k 
BL Po 3 = [E- 


=e 1 1 7; a —, c2 ; GS 
p- = = /—— | Foo — FF.) dy dy [19] 
16 7202 Aen Ti eae ieee 


jae 1 1 TRULY I"m 7 Pz 
Ba 16 72p,c° J r2 rs ee 
of: 


" T ret) 
= 6G, Se ooves os 2) Ge We | aap ain (ISI) 
i 2 ti \ lm Ue ij “Im | %Y %Y 


F dy dy” [18] 


and 


r? r4 ts i i 
of r.T.t ri? 

RSI gee gs SHES T dy d [21] 
74 74 Pe ij ~lm ij y vy 


Integrating the intensities in Equations [18] and [20] over the surface of a sphere gives 
the sound power radiated. The induction near-field terms in Equation [20] average out to zero 


over the sphere and thus do not represent radiated sound. The sound power radiated in the two 


cases is 
1 4n =, 3 
Be cogeallan FR dy dy [29] 
167*poc> 3 
and 
P : { a if, oF i) ay ay’ [23] 
hs era ra Saar, 90 95 aR 4 oe n0 y Yy 
2 5 15 BBW) HSER, 
16 7*po¢ 


(Equation 21, Reference 1) 


To examine the nature of the sound field due to a small volume of turbulence, it is 
helpful to change to spherical coordinates centered in the small turbulent region. For con- 
venience, 6 will be measured from the z,-axis and ¢ will be measured from the z-axis in the 


z,a,-plane. If r;=«a,, the following relations obtain for the components in the directions of 


increasing 7, 0, and ¢, respectively: 


7 0 @ 
Fa Oy oy 
—=sin@cos¢ — =cos¢ cos@ —=-singd 
Tr 0 p 
"9 9, oP) 
—s= sin@sin d —=sin¢g cos 0 —=cos¢ [24] 
r 0 ce) 
"3 05 i $3 
— = cos@ —=-sin@ —=0 
r 0 7) 


where 0 = |6| and ¢ =|¢|. 
These expressions can be used to show that, for the dipole case, the intensity is 


purely radial; namely, 


1 1 valle <8 : 
all cee dg dy [25] 
16 77p,.c° 2 72 


For the quadrupole case, however, the intensity has components in the direction of increasing 


7, 0, and ¢ as follows: 


He eye Bes cele T., 7,’ dy dy’ a 
1617p, {j OMe ia oe if 
and 
eal SDN or coin San ie T,’ dy dy’ 
16 7p,c [ p22 rd - 


In order to facilitate their use in actual cases, Equations [19], [21], [25], and [26] have been 
expanded and are given in Appendix A. Some physical insight is obtained if.the equations 
are examined in their vector form. For this reason, the equations for the mean-square sound 


pressure and the intensity are given in vector form in Appendix B. 


SPECIAL CASES 


The general sound field is too complicated to visualize readily, so in order to gain 
insight intothe general sound field a few special cases will nowbe examined. Four cases 
will be examined in some detail: (1) single dipole; (2) single lateral quadrupole; (3) single 
longitudinal quadrupole; and (4) quadrupole radiation from isotropic turbulence. Here the 
word ‘‘single’’ signifies a distributed multipole of one type and one orientation rather than a 
point multipole. A discussion of the changes in the spectra as determined by pressure meas- 
urements at various distances from these sources will be found on page 17 of this report. 

The sound field of a single dipole oriented along the z-axis will be discussed first. 
Substituting 7 = 7 = 3 into Equations [19] and [25] gives 


p? —afs — cos 20(Ris SRR) over [27] 
162 


and 
1 Le a congemag th, a? 
ee | — cos?6 FF,’ dy dy [28] 
16 7° poe 


The intensity is purely radial and has no induction near-field component. 


For a narrow band of frequencies around the angular frequency » in the spectrum of F, 


or in the sound radiated by the dipole, Equation [27] may be written in the form 


aD 1 1 iekroorarr ce? 
D -— | — cos Ou Er 1+ dy dy’ = (29) 
16 72c? r2 202 


Tr @ 


This indicates that there is an induction near-field term which changes the usual 1/r? varia- 
tion of the mean-square sound pressure when the square of the quantity, wavelength divided 
by 277, is significant in comparison with one. The sound intensity shows no change 
in variation with r upon entering the induction near field. In any plane containing the axis 
of the dipole, the intensity pattern is the familiar ‘‘figure eight”’ of cos? 6. 

The sound field of a lateral quadrupole oriented along the z,- and z,-axes will now be 


examined. Equation [21] gives the mean-square sound pressure for this case to be 


Fee \ 4 sin? cos? de0s? ( Fi Boy g fF Ts 
eae is N3isali3' 2B ibe} 34S} 
pic i ale 
+9 2 iF 7 dy dy’ [30] 


Similarly, from Equation [26], the sound intensity for this case is 


4 dei Een eee 
= rare sin? @cos’ cos? d T3713 ty ty [31] 
167° poe r 


The intensity is purely radial and has no induction near-field components. 


For a narrow band of frequencies around the angular frequency in the spectrum of 


T, or the sound radiated by the lateral quadrupole, Equation [30] may be written as 


2 4 
+9 


Tr @ rT @ 


== 4 tL ocr 
2 DP 2 2 
= any as i 0 0 1 
Dp \f 2 13213 Sin® 9cos* Ocos” } ( +3 


dy dy’ (32) 
16774 


The polar patterns of all terms in the mean-square sound pressure and the intensity 
are the same for the lateral quadrupole. In the z, z,-plane the pattern has the “‘clover leaf”’ 
shape shown in Figure 2. On the surface of any cone of constant 0 the pattern is a ‘‘figure- 
eight’? shape. The three-dimensional pattern may be described as four *‘beaver tails.”’ 

The sound field of a longitudinal quadrupole oriented along the z,-axis will be exam- 
ined next. The mean-square sound pressure, by substituting ¢ = 7 = 1 = m = 3 in Equation [21], 


is 


10 


o° 


6 
4sin®@ cos®@ 


Figure 2 — Directional Pattern of Intensity 
and All Components of Mean-Square 
Sound Pressure in Plane of the 
Two Axes of a Lateral 


90° 
Quadrupole 
180° 

= 1 iL Aico oe? 4 2 ear 
p= ae oi cos 6T,, ours (3 cos" 6 - 4cos*6+1) T;, T3, 

oc See 

+—(9c0s* 6 - 6 cos? 6 +1) T3373, | ¢y ay" [33] 
A 


Similarly, from Equation [26], the sound intensity for this case is 
1 1 A Seon Ge bs rer men get aia 
[= Sara he cos*@ 733 Ty, + (1 - 800s? 6) T3, 753 | 47 a 
16 7° p oe r r 


1 1 c2 : oo lee Le SBA 
Fearneln (oe Tilt) ee o 
0 
and 


Ig = 9 


For a narrow band of frequencies around the angular frequency » in the spectrum of 


T,3 or the sound radiated by this longitudinal quadrupole, these equations may be written: 


11 


c2 


(3cos* 6 — 4cos2 6 + 1) 


4 
No Fi (c0549 ~ cos24 +0] dy dy’ (35] 
rw 
Li Al Uisenrremms 4 c? 2 = i 
I= alae aa lng COSY tt (1 - 3cos* 6) | dy dy 
16 77p,c° re rw? “ 
eee <(-2sin@ cos 0)| dy ay” [36] 
Ree a ey ee —2sin 6 cos y ay 
16 7*p,c° 7? 1? Gy? 
and 
ly S 0 


If three equal, mutually orthogonal, longitudinal quadrupoles are combined, the mean- 


square sound pressure is 


ao 1 iL 1 20 20"? mee faba 
PAs merase T.,T,, a a [37] 
7 Cc 


As was pointed out by Lighthill,’ these represent the equivalent applied fluctuating pressures 
in the turbulent region. The induction near-field terms in the intensity also disappear. 


Figure 3 gives the polar patterns for the three terms in Equation [35] for all planes con- 


1 
taining the z-axis. The two induction near-field terms have zeros when @ = cos | — = 55 deg. 
3 


The middle term has a negative value for cos” 6 > 1/3 but at no distance or angle from the 
quadrupole does this make py negative. It may be seen that the induction near-field terms 
always predominate at 0 = 90 deg and the far-field term always predominates at @ = 55 deg. 

The polar patterns of the intensity due to a single longitudinal quadrupole are given 
in Figure 4, where the arrows indicate magnitude and direction. The far-field intensity is 
purely radial. The term for the induction near-field intensity, however, indicates a flow of 
sound energy out perpendicular to the axis and in again along the axis of the quadrupole, 
giving no additional net outflow of energy or power radiated. The induction near-field inten- 
sity is perpendicular to the surface of the cone generated by the 6 = 55-deg line. 

The sound field of a region of ¢sotropic turbulence shows no angular variation, but the 
induction near-field effects are still present. The assumption of isotropy permits the elimina- 


tion of many of the terms in such quantities as ih iN . To affect the simplification, it-will 
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Figure 3 — Directional Patterns of Three Components of Mean-Square 
Sound Pressure in All Planes Containing the Axis 
of a Longitudinal Quadrupole 
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Figure 4 — The Induction Near- and Far-Field Components of the Intensity 
Field in All Planes Containing the Axis of a Longitudinal Quadrupole 


Arrows indicate magnitude and direction. 
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be assumed that: 


T;; = Po%; (y, t — r/c) u; (7, t — 2/¢) = pot, u, [38] 
u,u,=0 if i= (isotropic turbulence)** [39] 
uu; uu, =U; uuu, + U, Uy u; Un + uu,” uu [40] 


(normal probability distribution)? 


and 


u,, =-U,U,, u,u,;=0 ete. [41] 


(stationary random functions of time) 


Also, let & - 7/r = u be the component of the turbulent velocity fluctuation w@ in the 7-direction 
at point 7 and at retarded time (¢ - r/c). The velocity u is also, on the average, equal to the 
component of w in any other direction. 


By using these equations, Equation [21] reduces to 


2 
“Ty. Po 4 —— Tom c? Peery cn ie o Saar, 
p = ff Af farar. ay} +2— wu’ wu’ +3 —(uu’)*fdy dy’ =: [42] 
16 77c* 7? 7? is 


r 


Similarly, the intensity from Equation [26] reduces to 


2 

Po 4 {|__ __, SaanD 5 

I,- ——|| —|uu’ ti’ + 3 (au’)*| dy dy [43] 
1677p ,crry 


The 6 and ¢ components of the intensity are zero. Equation [43] corresponds to the result 
of Proudman, Equations 4.5 and 5.1,3 with the ‘‘decay’’ terms neglected. 

As was pointed out by Lighthill,! the spectrum of the sound radiated is different from 
that of the velocity fluctuations because the radiated sound pressure depends on the square 
of the velocity fluctuations. Therefore, the sound pressure will include sum and difference 
tones of the velocity frequencies and, in general, will have a flatter spectrum. If the veloc- 
ity fluctuations are confined to a narrow band of frequencies around the angular frequency w* 
one might approximate the effect of sum and difference frequencies by assuming that the angu- 
frequency o of the radiated sound is approximately equal to 2w*. The mean-square sound 


pressure and intensity under this assumption then reduce to 
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=> Po 1 Z 4 
p= \{= w* (uu’)? (1 +9 + 12 ay dy’ _—s« ([44] 
16 72c* r2 r2@2 ro 
or 
3 Po 16 2 4 
2= le (win’)? (1 +2 + 12 ag dy’ [45] 
1672c* 7? ra rte 
and 
Po Gis 
I, = a eles (ui)? dy dg’ [46] 
16 mpc 72 


As one would expect, there is no angular variation in the sound field, but the mean-square 
sound pressure has two induction near-field terms. This sound field turns out to be equivalent 
to the individual contributions of a longitudinal quadrupole in each of the coordinate direc- 
tions and two lateral quadrupoles in each of the pairs of coordinate directions. 

The modifications in the equations required when the flow is analyzed with respect to 
a moving frame have been discussed by Lighthill? and Phillips‘ for the far-field radiation. 
Similar modifications may be made for the general sound field. Because of the complexity of 
the equations, however, this has been done only for the instantaneous and mean-square sound 
pressures. These will be found in Appendix C. 

The utility of the foregoing theoretical results depends on the ability of observers to 
estimate certain, up to now, unmeasured quantities. These estimates may be made more realis- 
tic by using a number of useful concepts and by keeping in mind the special properties of the 
near-sound field. These concepts will be discussed in the next section for a few practical 
applications. 


DISCUSSION 
CORRELATION VOLUME AND CORRELATION AREA 


For purposes of estimating some of the integrals introduced previously, the concepts 
of correlation volume and correlation area are useful.1’?’4*> The correlation volume and corre- 
lation area as introduced by Lighthill refer to an average eddy volume and area inside of which 
the fluctuating quantities are considered correlated and outside of which the fluctuating quan- 
tities are considered uncorrelated. The correlation volume and correlation area may be given, 
respectively, as 


TT,” 
L 
V ae dy’ [47] 
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and SS 


ae 
S -| wiz) [48] 


These expressions also hold for the time derivatives of the quantities under the integrals. 
Since eddies of different sizes have different frequency scales, V and S defined for a narrow 
band of frequencies would be expected to be functions of the frequency. One would expect 

V and S to be smaller for the higher frequencies; that is, that V would be approximately pro- 
portional to the inverse cube of the frequency and S approximately proportional to the inverse 
square of the frequency. This exact proportionality, of course, neglects the effects of sum- 
and-difference frequencies in the velocity fluctuations and ignores the fact that time fluctua- 
tions are not synonymous with the space variations that are usually measured. The frequency 
dependence of V and S partially offsets the higher efficiencies of quadrupole and dipole radia- 
tion at the higher frequencies. As a result, the frequency of the peak in the sound spectrum 
should be about the same as the frequency of the peak in the spectrum of the fluctuating 
velocities or fluctuating forces. The correlation volume and correlation area might also de- 
pend on the particular velocity components and force components used in the integrals. In 
spite of these uncertainties, the effect of integrating over dy’ may be approximated by replac- 
ing dy’ by V or S and removing the primes from the integrals that give the mean-square sound 


pressure and intensity. 


MEAN-SQUARE SOUND PRESSURE IN TERMS 
OF SOUND POWER RADIATED 


The mean-square sound pressure at any point in the sound field outside the geometric 
near field of like-oriented dipoles or quadrupoles may also be given in terms of the sound 
power P radiated by the multipoles. For the four special cases discussed previously (namely, 
dipoles, lateral quadrupoles, longitudinal quadrupoles, and a region of isotropic turbulence), 


the mean-square sound pressure is given, respectively, by 


Lo Py Gl? 2 
2 2 
Da = (3 cos~ @) ( + [49] 
Agr2 reas 
1) MOa@le 2 of 
p? = (15 sin? 6 cos? 6 cos” ¢) (108 +9 ) [50] 
4nr? ro Ge 


pocP 


2 4 
= c 4 1 c 6 1 
p? = (5 cos* 6) + - + )s (s - + | [51] 
Aare rq? cos?@ cos*é ret cos?6@ cos*é 
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The quantity in each case gives the mean-square sound pressure at all points in the 


4ar 
sound field of an equivalent simple source. The terms that depend on 6 and ¢ indicate the 


variations of p” from the average over all directions. The last bracket in each equation gives . 


the induction near-field correction to the radiated mean-square sound pressure. 


SLOPE OF FREQUENCY SPECTRUM 


It is apparent that the effect of the induction near-field sound pressure is to enhance 
the low-frequency components in the sound spectrum as it appears in the far field. Therefore, 
the sound pressure spectrum determined from data on the induction near field will differ from 
the spectrum determined from data on the far field. The effect of the induction near field is 
to decrease the slope of the spectrum by 6 db/octave in the case of dipoles and 12 db/octave 
in the case of quadrupoles. For a narrow band of frequencies the approximate far-field sound 
pressure level may be obtained from the induction near-field sound pressure level by sub- 


tracting the quantity 


2 
72 rt e4 


ce? ct 
10 log ( + a, +a } [53] 


in each case using a mean frequency of the band and an average distance 7. The quantities 
a, and a, are to be taken from the last bracket term in Equations [49], [50], [51], and [52] for 
the four cases discussed. The application of these equations to the study of the sound from 


turbulent jets and turbulent boundary layers will now be discussed. 


SPECIAL CASES 


The case of a turbulent jet is of special interest because of its widespread application. 
The sound produced by a jet is usually considered to consist of two parts: (1) from a cylin- 
drical region of strong shear near the exit of the jet; and (2) fromthe core of the jet and a 
larger mixing region downstream. The sound from the shear layer is predominantly from lateral 
quadrupoles with axes along the jet and perpendicular to the shear layer of the jet. Thus, if 
6 is measured from the axis of a circular jet and if the observation point is at least several 
jet diameters from the jet, the term cos“ ¢ may be replaced by its average, namely 4, in the 
equations pertaining to lateral quadrupoles. As a first approximation, it may be assumed that 
the sound from the core and the downstream part of the turbulent region of the jet is given by 


the equations for the quadrupole radiation from isotropic turbulence. 


17 


To obtain reasonable corrections for data taken in the near-sound field of a jet, an 
assumed distribution of quadrupoles along the jet is integrated with respect to varying dis- 
tance and angle from the measuring point. This procedure automatically includes the effects 
of the geometric near field, the induction near field, and the angular dependence of the quad- 
rupoles. To a first approximation an average distance and angle may be used. For high con- 
vection Mach numbers it may be necessary to use the equations in Appendix C. 

For the typical sound spectrum of the far-field sound of an air jet, given by the solid 
curve in Figure 5, the approximate spectrum expected from data taken at a distance of 1 ft 
from the jet is indicated by the dashed curve in the figure. The frequency corresponding to 
a 1-ft wavelength in air is about 1100 cps. The eventual change in slope at the lower fre- 
quencies is 12 db/octave, but the change in slope in the transition region between the near 
field and the far field is somewhat indeterminate because of the assumptions made as to the 
average distance 7 and the types of quadrupoles. These estimates can be made sufficiently 
accurate, however, toenable one to make reasonable corrections to data taken inthe near- 
sound field. 

The case of sound from a turbulent boundary layer is also of interest. Aspects of this 
problem have beendiscussed by Lighthill,?* Powell, > Phillips, * and Curle.” Curle has 
shown, in general, how dipole radiation results when a solid body or boundary is present in 
the turbulent field. For curved boundaries dipoles with axes perpendicular to the boundary 
arise from the fluctuating pressures on the boundary, and dipoles with axes parallel to the 
boundary arise fromthe fluctuating shear forces. The dipoles from the fluctuating forces 
exerted by a plane boundary on a fluid have been shown by Phillips to have their axes parallel 
to the boundary. If these dipoles are assumed to be of approximately equal strength in the 


two coordinate directions in the surface and if, for convenience, the origin of coordinates is 
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Figure 5 — Typical Air Jet Noise Spectrum and Approximate Changes that Occur 
if Data Are Taken at an Average Distance of 1 Foot from the Jet 
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located at the measuring point, the contributions from a given area of uniform boundary layer 
may be easily found by integration. For instance, by using Phillips’ tentative estimate for 


an z-Reynolds number of approximately 108 (sound power radiated per unit area equals approxi- 


2 
ce 
mately 4(1077) Po oP ey and by using Equation [49] with P(8cos? 6) (: + \e 


rT @ 


replaced by 
3 3 , , 
6, 8(10°’)pyU'MS (; sin?) cos? 6 hoy eeaia. ae 
2 
| 1+ —— —— 


; h? 
Oy 


de’ 
cos? 6’ 


the broad-band mean-square sound pressure a perpendicular distance / from the center of a 


circular ring area between 0’°= 6, and 0°= 0, is 


2 


cos 6, a 


, 


5 1 
p? = 6(10-’)p,-U*M [-( ) A (cos? 6, - cos? 6,) + (sin* 0, - sin* 00] [54] 


cos 6, 40h? 
where M, = U/c, 
U_ is the free-stream velocity, and 


6 is the angle with the perpendicular to the plane boundary. 


Phillips has shown further, however, that the dipole contribution is negligible if the boundary 
layer is homogeneous in layers parallel to the boundary. Thus, for this condition the predom- 
inant sound would be from quadrupoles in the surface and in the turbulent boundary layer. A 
procedure similar to the above might be used to estimate the contributions of these quadrupoles. 
The quadrupole radiation from the turbulence proper might be treated approximately as contri- 
butions from isotropic turbulence outside the shear layer and from lateral quadrupoles in the 
shear layer just outside the laminar sublayer. Further, the effects of aconvection velocity 
may be included by using the equations in Appendix C. 

The above discussion of the noise from a turbulent boundary layer holds only for points 
in the fluid outside rigid bodies. In the case of nonrigid bodies, motion of the boundary, in- 
duced by the fluctuating pressure field on it, will create a system of simple sound sources 
over the boundary and will permit propagation of dipole and quadrupole sound through the 
boundary. A thin pliant boundary with turbulent fluid motion along one side and no net fluid 
motion along the other side, for instance, might radiate simple source, dipole, and quadrupole 
sound both intothe moving fluid and intothe motionless fluid, withthe multipole sound radiated 
into the motionless fluid modified by the transmission properties of the boundary. The simple 
sources do not have an induction near field but they can have a geometric near field. The 
transmission and reflection properties of a nonrigid boundary in the induction near field of 
dipoles and quadrupoles are, of course, not as simple as for a plane sound wave. Because 


of their greater efficiency, the simple sources should predominate at low Mach numbers for 
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many actual cases of turbulence on nonrigid boundaries. Any induction near-field correction 
thus will depend significantly on the relative contributions of the sources, dipoles, and quad- 
rupoles in the boundary layer. To facilitate the determination of these relative contributions, 


the effects of the Mach number and other factors should be considered. 


DIMENSIONLESS RATIOS 


A little dimensional analysis facilitates the understanding of some of the properties 
of the induction near-sound field. By using a typical velocity U and a typical linear dimension 
L in the flow, as in Reference 1, and by ignoring the numerical factors, the expressions for 
the ratio of the mean-square sound pressure to the square of the typical Bernoulli pressure 
for the source, dipole, and quadrupole, respectively, may be written 


2 2 
D L 
a 5 
ji V2 5 [55] 
Bo 
p? L2 14 
a a me [56] 
1 2 2 72 0 r* 
— pU 
and 
yD) 2 4 6 
IL L iG 
iA ETE A ete a ease Ng 2 es [57] 
1 9\ 2 7? r4 7? 
Ope 


These expressions indicate that very near the turbulent region, where L/r is of the order of 
one or more, the root-mean-square sound pressure is of the order of the typical Bernoulli 
pressure in the flow itself and independent of the Mach number at low Mach numbers. This 
also indicates that very near the turbulent region the density changes are of the order of the 
density changes in the flow itself. 

It is of interest also to know the ratio of the absolute value of the fluctuating Bernoulli 
pressure to the absolute value of the sound pressure at any point in the field of sources, di- 
poles, and quadrupoles. In the far field one can show that this ratio is always equal to the 
ratio of the particle velocity to twice the sound velocity or 

i 2 
9g Po? | ISK ay| 
ee ee [58] 
p 87 Po re? 


for a source, 
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ee Beh 
9 Po | fF, dy| cos 0 


[59] 
P 8 7p ytc> 
for a dipole, 
i 2 = 
Or | fT,, ¢7| sin 0 cos 0 cos¢ 
RS NRT TSS Ee [60] 
y 8 TP ote 
for a lateral quadrupole, and 
1 2 aA 
%@ Po" ities dy| cos? 6 
a Ni 161] 


4 
P 8mpyre 
for a longitudinal quadrupole, where K = K (y,¢-7/c) is the ‘‘mass rate density’’ or acoustic 
source strength per unit volume of the sources. 


In the induction near field the ratio is equal to the ratio of the particle dis- 


1 2 
reanolg /P 


placement to twice the radial distance from the multipole, with an additional angular depend- 
ente for the higher-order multipoles. In the induction near field of the four cases above, the 


ratios are, respectively: 


1 
fs p v2 = 
9 Po lf K dy| 


= —______——_. [62] 
P 87p,o" re 
1 2 4 
Bw If F, dy] (1+ 3cos? 6) 
$$ |= —<<$<___. ———_—_ [63] 
p Brp yor |cos 6| 
1 
2 se 2 : 
9 Po” |fT,, dy| 319 sin? 6 cos? 6 cos? ¢ + cos* (cos? 6 - sin? 6) + cos*@sin? 4] 
Dp 87p,o* ° |sin 6 cos 6 cos d| 
[64] 
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1 
2 oA 2 
% Po” | fT3, @| 9[(1 - 3cos* 0) + 4cos?@ sin? 6] 


Ny aU TT eT ea [65] 


p Bmp ,o'r |3cos? 6-1| 


In general, the sound pressure from randomly oriented multipoles predominates over the fluctu- 
ating Bernoulli pressure in the induction near field as long as the magnitudes of the particle 
motions induced by the multipoles are small compared to the distance from the multipoles. 

The dependence of the ratio of the sound power radiated to the mechanical power 
dissipated (the so-called ‘‘efficiency of conversion of mechanical power to sound power’’) 
on the Mach number and other factors may also be used to estimate the relative contributions 
of sources, dipoles, and quadrupoles. Except for secondary effects, this efficiency is propor- 
tional to M, for sources, Mo for dipoles, and Mo for quadrupoles. 

For jets the efficiency of conversion of mechanical power to sound power has been de- 
fined in terms of the total mechanical power of the jet. The efficiency of conversion to sound 
for aboundary layer, however, is more conveniently expressed in terms of the time rate of 
energy dissipation by the boundary layer; namely, the drag times the free-stream velocity. In 
terms of the drag coefficient Ces the efficiency of conversion to sound for the dipole and quad- 


rupole radiation from a boundary layer may be written, respectively: + 


acoustic power per unit area 
0g = Mg (Bs Mo) [66] 
ie 3 
5 Pol Cr 


acoustic power per unit area M590. (R [67] 
ur ay 1 3 ee, Iq (R, io) 
Baye 


Thus, in ranges of conditions in which each multipole source predominates, data may be con- 
solidated and 7, and ily evaluated. These efficiencies would be expected to show secondary 
Mach number effects and Reynolds number dependences through the functions g,(#, Mg) and 
Iq (FR, Mo), especially at low Reynolds numbers. 

For a uniformboundary layer, the rate of energy dissipation bythe layer should be a 
more realistic standard thanthe total power of the jet as used to compare jet noise data. 
Thus, the efficiency 1g would not be expected to be the same as the acoustic efficiency as 
defined for a jet. The presence of high viscous dissipation of energy very near the boundary 
would tend to make 1g lower than jet efficiencies. But the fact that much of the energy dissi- 
pated in a jet is in a region with lower Mach number than that used to define the jet sound 
conversion efficient, would tend to make Ng higher than jet efficiencies. The separation of 


the contributions of the fluctuating forces and the fluctuating stresses and the evaluation of 
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Nd and 1, may be facilitated by using the special properties of the induction near-sound field 


of dipoles and quadrupoles. 


SUMMARY 


In summary, the near-field properties of the general sound field of dipoles, lateral quad- 
rupoles, and longitudinal quadrupoles, as well as isotropic turbulence, have been discussed 
for various assumed conditions. The general expressions for the rms sound pressure are ex- 
pected tobe the most useful in practical applications. These applications would include the 
analysis of internal noise and structural loads caused by jets, propellers, and the boundary 
layer of aircraft and ships; the analysis of surface and turbulence noise in the close vicinity 
of wind tunnels and water tunnels; and the analysis of surface and turbulence noise in the 
vicinity of the boundary layer of objects moving in air or water. To facilitate the use of the 


equations inthese applications, a number of useful concepts are also given. 
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APPENDIX A 


EXPANDED FORM OF EQUATIONS FOR MEAN-SQUARE SOUND 
PRESSURE AND INTENSITY 


Equations [19], [25], [21], and [26], respectively, may be expanded by letting ¢, j, 2, 
and m take on the integral values of 1, 2, and 3 in all combinations and summing the terms. 


These equations then become 


see Be 
p= edie (44° fy ia dy dy’ [68] 
16 77c 72 
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eS eto 
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where 
A =sin@ cos¢ F, + sin@ sing F, + cos 6 F, [72] 


aD 2 5 my 0.9) 2 
B = sin“ 0 cos ¢ T,, + sin” 6 sin” ¢ T,, + cos 6T3, 
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DVO Dee rare [74] 
E = sin @ cos 0 cos? d T,, + sin @ cos @ sin? ¢ T,. — sin@ cos @ T,, [75] 
F =2sin@cos@sin¢cos¢ T,,-(1-2c0s70) cos ¢ T,3-(1-2cos76) singé T,, [76] 


G =-sin@ cos¢ sind T,,+ Sin9 cos¢ sing T,, [77] 
and 

H = sin (1-2sin? 4) T,, - cos@ sing T,, + cos@ cos¢ T,, [78] 
It may be seen that, in general, terms involving interaction between lateral and longitudinal 


quadrupoles, longitudinal quadrupoles in different directions, and lateral quadrupoles in dif- 


ferent directions, are present. 
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APPENDIX B 
VECTOR FORM OF EQUATIONS FOR MEAN-SQUARE 
SOUND PRESSURE AND INTENSITY 


By defining the fluctuating force vector F by the relation FF, = F . F, Equations [19] 


and [25] may be written in the vector form 
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By using Equation [38] and the relation u,u,; = @ - u, Equations [21], [20], and [26] may be 


written inthe form 
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APPENDIX C 
EFFECT OF CONVECTION OF TURBULENCE 


For the case of fluctuating forces being convected at vector Mach number M the trans- 
formations used by Lighthill! and Equation [1] lead to 


1 rF, Q1-M?)7,F, MAF. 
2 esd ll ees ee oo Smee nome (td [84] 
47) | o(r-M-7F)2 (r-M-F)3— (9 - MF)? 


where 7 = 7 + Mr and F, = F, (j, ¢- 7/c). By taking the time average of the square of the 
instantaneous sound pressure, the mean-square sound pressure becomes 
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where terms involving odd numbers of time derivatives have been neglected as before. 
If the fluctuating forces are being transported in the z-direction, |M| =M = M, and 
M.7/r=Mcos 6. Therefore, the mean-square sound pressure is 
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[86] 
For the special case of dipoles having axes in the 2,-direction and being transported at Mach 
number M in the z,-direction, the mean-square sound pressure is 
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Similarly,for the case of turbulence being convected at vector Mach number M, the in- 
stantaneous sound pressure, by using Equation [2], is 
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where now T;; = T;, (7, t-r/c). By again neglecting terms involving an odd number of time 
derivatives, the mean-square sound pressure becomes 
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If the turbulence is being transported in the 2,-direction, the mean-square sound pres- 


sure becomes tr 
3(1-M2)- 
ay 1 ij 1 Crh em Bae ec 3(1- rTP 
: 167204) 72(1-M cos 0)2 |r4 (1-M cos 6)* 7” = 72 (1-M cos 6)? } | r4(1-M cos 6)4 
4(1-M) 5,1) P, A eae 12M(1-M?) 2, 1.7, = 16 M?r, 7, a 
+ OF Lo Dy iis (oer arererarieee nd Ua Bri ep cea ere ee hey Uc 
72(1—M cos 6)? Pee Al 73(1-M cos @) ‘ r*(1-Mcos 6)” i 
2 dl 2,4 
4M rr = 8Mr, — of 9(1-M*) retry 
= T,, 1. Aes Ty a ee 
r2(1—M cos 6)” r(1—¥ cos 6) r*(1—M cos 6)* 4 (1—Mcos 6)* 
2,3 on 
6(1-M*) 6..7,7 386(1-M*) r. 7.37 
ij lim 2 2 5 a7 eb ——_ 
Se ae le) oF; 1m Tj 1B Ver sora a Ta , T)3 
r?(1-M cos 6)? 73(1-M cos 6)3 
2 2? 2 2)? 22 
36 M2(1-M?) rr, 12M7(1-M) 7, 7, 12M(1-M?). r, 
+ ———- T., T,, + ——__ T_T, ¥ + ——_ 7. T’ 
3 413 33 3 *1l 
r*(1-M cos 6)2 ‘ r?(1-M cos 6)? i r(1—M cos 6) : 
24. M3 (1—M?) 7, slates Eee 
SS TES — IE (LSI) Tp TES 2 OI RS TBS (CI I? [90] 


r(1-M cos@) 


For the case of lateral quadrupoles having axes in the #,- and z,-directions and being trans- 


ported in the z,-direction at Mach number M, the mean-square sound pressure is given by 


ar 1 “(¢ 4 sin? 6@cos?¢ Sapa 
2 2 , 
\| a GOS Hn US 


Pp = 
16 72c4 r7(1-M cos 6)® 


2 Petia mere 

ce? |38(1-M?) cos?6@ 6M(1-M?) cosé Eli cece 

ee ti | 7 Ts 
2 2 

r (1-M cos 6) 1-M cos@ 


4 3 2 
c* 19(1-M?) cos?@ 18M(1-M2) cosé : 9M? (1-M?) 
r4 (1—M cos 6)* (1-M cos 6)3 1-M cos 6)2 


T13T;3 dn dn” 
[91] 


Beyond several diameters from a circular jet, cos? may be replaced by its average value; 
namely 44, as before. 
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For the case of longitudinal quadrupoles having their axes in the z,-direction and 


being transported in the 2,-direction at Mach number M, the mean-square sound pressure is 


given by 

a 1 1 Sense 

pe = \\ ‘Titi yh Mam Me 
16 72c* 72(1-M cos 6)4 

2 
92 |3(1-M?) cos*@ 12M(1-M”) cos* 6 ps 
aay ——E_ (a) eoseioOlieostarl Ta ie 

r (1-M cos 6) 1-M cos @ 


4 4 3 2 

©’ | 9(1-M2) cos*@ 36M(1-M2) cos?6 6(1-9M7) (1-M?) cos? 6 
P|) ee Se eee 

ro (1-M cos 6)* (1-M cos 6)3 (1-M cos 6)? 


12(1-3M?) (1-M?) M cos 6 3 = 
+ (1-3M Tina (aniann 92 
1-Mcosé 33/538 ee ee 


It may be seen that introducing a convection velocity somewhat complicates the ex- 
pressions for the general sound field. The intensity could be obtained for the case of convec- 


tion of dipoles and quadrupoles, but the complexity would be great indeed. 
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